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1.  Introduction 

Active  control  of  sound  propagation  in  ducts  plays  an  important 
role  in  many  technical  applications,  e.g.  the  attenuation  of  noise 
in  ventilating  systems  or  in  the  exhaust  systems  of  combustion 
engines.  Today  many  control  strategies  exist  for  the  simplest 
case,  which  is  the  noise  control  in  the  low  frequency  range  in 
ducts  the  lateral  dimensions  of  which  are  small  compared  to  the 
wavelength  of  sound.  Up  to  the  cut-on  frequency  of  the  first 
higher  order  mode  a suppression  of  the  propagating  plane  sound 
waves  is  possible  in  a broad  frequency  range. 

The  situation  is  more  difficult  when  the  lateral  duct  dimensions 
become  large  compared  to  the  wavelength,  as  for  example  in  the 
intake  section  of  an  aircraft  engine.  With  increasing  rotor  speed 
an  increasing  number  of  acoustical  modes  becomes  present.  This 
means,  the  system  comprises  an  increasing  number  of  degrees  of 
freedom.  As  a consequence,  control  of  the  sound  field  is  practi- 
cable just  in  a narrow  frequency  band  or  even  limited  to  discrete 
frequencies. 

In  figure  1 the  frequency  spectrum  of  the  sound  field  in  the  duct 


section  of  a turbo-machine  is  depicted.  Obviously  the  fan  pro- 
duces tones  of  high  intensity  compared  to  the  broad  band  noise 
level.  Attenuation  of  these  tones  by  means  of  passive  absorbers 
often  is  not  sufficient.  Therefore  the  active  suppression  of  the 
tones  represents  a reasonable  improvement. 

A common  approach  for  active  noise  control  consists  in  the  use 
of  the  well  known  adaptive  LMS-algorithm.  This  strategy  ig- 
nores to  a certain  extent  the  properties  of  the  physical  system  to 
be  controlled.  Furthermore,  with  increasing  complexity'  of  the 
modal  sound  field  the  numerical  costs  of  the  adaptive  system 
increase  fast  and  the  problem  of  numerical  instabilities  occurs. 
Thus  a physical  approach  is  reasonable,  which  aims  at  minimi- 
sation of  only  the  dominant  modes  of  the  undesired  tone.  This 
results  in  a reduction  of  the  control  effort.  Furthermore  the 
stability  of  the  complete  modal  control  process  is  improved  due 
to  its  separation  into  a number  of  independent  subprocesses  for 
individual  modes.  Besides,  as  another  control  strategy,  the  di- 
rectivity of  the  radiated  tones  in  the  far  field  outside  the  duct  can 
be  influenced  through  a specific  control  of  the  modal  structure. 


Figure  1:  Sound  pressure  spectrum  detected  by  a micro- 
phone in  the  duct  section  of  a turbo-machine  (bpf=blade 
passing  frequency). 


In  this  regard  a azimuthal  mode  control  strategy  was  developed 
and  experimentally  proved  by  NASA  Lewis  Research  Centre  in 
1996  (Fla  and  Hu  [1]).  By  use  of  a single  ring  with  pressure 
sensors  and  one  ring  with  actuators  they  were  able  to  reduce  one 
dominant  azimuthal  mode  at  the  axial  ring  location  of  the  sen- 
sors. 

As  the  amplitude  of  the  azimuthal  mode  is  a function  of  the 
axial  and  the  radial  position  in  the  duct,  control  at  a single  axial 
ring  position  docs  not  ensure  minimisation  of  that  azimuthal 
mode  at  other  duct  positions,  in  particular  if  the  sound  wave  has 
a high  radial  content. 

A direct  radial  mode  control  strategy  at  present  seems  to  be 
impossible  in  practice,  because  variations  of  temperature,  fre- 
quency and  flow  velocity  in  the  duct  strongly  influence  the 
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accuracy  of  the  online  radial  mode  analysis.  Besides,  a large 
number  of  variables  must  be  considered  by  the  control  system, 
leading  to  instabilities  in  the  minimisation  process. 

Therefore,  an  improved  control  strategy  has  been  proposed  by 
Zhang  et  al:  the  indirect  radial  mode  control  (IRMC),  which 
enables  effective  control  of  the  modal  sound  field  without  di- 
rectly determining  the  radial  amplitudes.  The  method  is  based 
on  the  well-defined  relationship  between  the  azimuthal  ampli- 
tudes at  different  axial  positions,  which  is  given  implicitly  by  its 
composition  of  a fixed  number  of  radial  components.  If  one 
succeeds  in  reproducing  the  azimuthal  amplitude  and  thus  the 
correct  linear  superposition  of  radial  components  at  different 
axial  positions,  suppression  can  be  achieved  not  only  at  this 
distinct  locations  but  also  in  the  whole  duct. 

In  the  following  chapters  it  will  be  shown,  that  the  performance 
of  the  modal  control  strategy  strongly  depends  on  the  proper 
choice  of  the  number  as  well  as  the  locations  of  the  sensors  and 
actuators. 

2.  The  sound  field  description  in  circular  and  annular 
ducts 

The  propagation  of  sound  pressure  waves  in  cylindrical  or  an- 
nular ducts  in  the  presence  of  unifonn  mean  flow  can  be  de- 
scribed by  a linear  acoustic  wave  equation,  see  e.g.  [2],  A solu- 
tion of  the  wave  equation  in  cylindrical  coordinates  ( x,r,<p ) is 
given  by  a mathematical  decomposition  into  azimuthal  and 
radial  components.  In  the  frequency  domain  it  can  be  written  as 
(see  [3]): 

p(x,r,(p,co)=  ’^Am(x,r,(0)-e""'p  (1) 

Am  ( x , r.  <0)=  Y /„„  (Om„  —)[a„  (CD)e,<;"'  + A~n  ((0)e‘t""x } (2) 

ti  n, 

where  A,r  (x,r,a)  is  the  amplitude  of  the  azimuthal  mode  of 
order  m.  Am(x,r,0) ) expands  into  a series  of  radial  modes  with 
amplitudes  A^(a>)  and  Arm(a>) . “+”  denotes  wave  propaga- 
tion in  flow  direction,  against  flow  dirction.  k*r{a>)  and 

K,m  (®)  are  the  respective  axial  wavenumbers,  r,  is  the  outer 
duct  radius,  and  fmn  denotes  the  cylinder  functions  which  de- 
scribe the  radial  distribution  of  the  modal  sound  pressure  field  in 
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Figure  2:  Longitudinal  section  of  the  duct. 


the  circular  or  the  annular  duct  (a  detailed  description  can  be 
found  in  [3]). 

3.  The  Indirect  Radial  Mode  Control  (IRMC) 

In  figure  2 the  intake  of  a turbo  machine  is  depicted  with  a 
sensor  and  an  actuator  array  mounted  flush  with  the  inner  duct 
wall.  The  sensors  and  actuators  are  arranged  in  several  axial 
rings.  In  each  ring  the  sensors  resp.  actuators  are  spaced  evenly 
in  the  azimuthal  direction  (p  . The  number  of  sensors  and  ac- 
tuators per  ring  must  be  chosen  according  to  the  Nyquist  theo- 
rem. The  arrangement  of  the  sensors  allows  a azimuthal  mode 
analysis  for  each  ring  as  was  demonstrated  in  [41.  It  is  a prereq- 
uisite for  the  IRMC-method  processed  in  the  azimuthal  mode 
domain.  Figure  3 gives  an  overview'  on  the  required  signal  proc- 
essing. 

With  the  azimuthal  modal  analysis  realised  as  a real  time  appli- 
cation (see  Bohning  et  al  [5])  the  modal  amplitude  Pm(xi,ra,0)\ 
for  each  azimuthal  mode  order  -m,  ....  +m  of  the  primary  sound 
wave  (produced  by  the  rotor)  can  be  determined  for  the  set  of 
sensor  rings  { }.  The  objective  of  the  JRMC  is  to  cancel  the 
azimuthal  modes  Pm(x,,ra,0})  by  a secondary  wave  of  adequate 
azimuthal  mode  amplitudes  Sm  {xnra  ,(0) , which  are  generated 
by  the  actuators.  Mathematically  this  is  equivalent  to  minimisa- 
tion of  the  sum  of  these  mode  amplitudes  over  all  axial  positions 
of  the  sensor  rings,  seperately  for  each  azimuthal  order  m: 

/»(®)  = XI  P’»(xj,r°’co')  + S'»(xi.r°’coi2  =Min.  (3) 

j=  I 

The  secondary  mode  amplitude  al  ihe  sensor  locations  { x,  } 
generated  by  the  actuators  is  given  by 

u 

S,„  (*, , ra,(0)  = £ (t,„  ),„.  (or)  • ( xh  , r„ , <u) , (4) 

h= 1 

where  q,„(x4,r  ,ft>)  represents  the  azimuthal  mode  amplitude 
of  the  actuator  ring  at  the  axial  position  xt  . (tm)hJ  is  the  coeffi- 
cient of  the  modal  transfer  matrix  7y  which  describes  the  mo- 
dal transfer  functions  between  the  axial  positions  { xh  } and 
{ x , }.  The  modal  transfer  matrix  has  to  be  determined  in  a 
separate  step  prior  to  application  of  ANC,  see  chapter  4. 

From  equation  (2)  it  can  be  deduced,  that  in  order  to  control  an 

primary  fan  noise 


Figure  3:  Schematic  of  the  signal  processing  for  IRMC. 
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azimuthal  mode  consisting  of  N radial  modes,  at  least  N sensor 
rings  are  necessary.  If  the  azimuthal  mode  amplitude 
P„(x.,r,i o)  at  these  N axial  positions  x.  is  reproduced  on  the 
duct  wall  (r  = ra ),  the  radial  structure  in  the  duct  is  also  repro- 
duced implicitly. 

For  simplicity  a vector  of  azimuthal  mode  amplitudes  of  mode 
order  m at  the  set  of  axial  positions  { x,  } resp.  { xh  } is  intro- 
duced: 

Pm  :=  [Pm(xiru,(0),Pm(xlra,(o),....,Pm(x,  ri:,co)j  , 

S_m  :=  [Sm(xl  r:l,(i)),Sm(x2  r.,co),....,Sm(Xj  ru,co)J  , 

lm  ■=  [!„(*, >!»(**/; >®)]'  • (5) 

In  equ.(5)  single  tones  with  constant  frequencies  m are  consid- 
ered. Since  only  wall-flush  mounted  sensors  and  actuators  are 
used,  r can  be  dropped  in  the  further  analysis.  Thus  equation 
(4)  can  be  rewritten  as 

Sm=ThJ-%  , (6) 

and  one  obtains  the  solution  of  equation  (3)  as  least  square  fit 
expression 

L, = _^=I  )'*  • tL ■ (I ! )'’•£»•  {1) 

Since  the  system  is  linear,  the  control  process  can  be  performed 
separately  for  each  mode  order  m of  a single  tone. 

The  final  actuator  signals  (in  the  frequency  domain)  can  be 
computed  according  to  equation  (1)  for  each  actuator  at  position 
(xk,rl,(pl) . In  order  to  generate  the  secondary  tone  with  the 
correct  phase  leading  to  the  cancellation  of  the  primary  field,  it 
is  locked  to  a trigger  signal,  which  is  produced  by  the  rotor. 

Theoretically  the  system  of  linear  equations  (7)  can  be  solved 
exactly.  In  practice  the  inversion  of  the  included  hermitian 
matrix  may  be  critical  numerically  and  lead  to  instabilities.  The 
problems  that  appear  in  connection  with  a matrix  inversion  will 
be  discussed  in  more  detail  in  chapter  5. 

4.  Simulation  of  the  IRMC 

Before  performing  IRMC  the  modal  transfer  matrix  has  to 
be  determined  by  an  azimuthal  mode  transfer  function  meas- 
urement for  each  azimuthal  mode  order  -m  ,...,  +m.  The  transfer 
matrices  include  the  acoustical  transfer  path  as  well  as  the  elec- 
tro-mechanical transfer  characteristics  of  the  actuators  and  the 
respective  signal  processing.  For  the  numerical  simulation  and 
analysis  of  the  IRMC  here  the  modal  transfer  matrix  Tv  is 
determined  theoretically. 


In  order  to  get  a more  realistic  azimuthal  modal  spectrum  as  an 
input  for  the  simulation,  a modal  noise  amplitude  appropriate  to 
a measurement  error  of  the  sensors  is  added  to  the  modal  am- 
plitudes in  Pm  (co) . The  magnitude  of  the  noise  amplitude  is 
related  to  the  maximum  modal  amplitude  by  considering  a given 
‘signal-to-noise-ratio’  a . 

In  figure  4 the  simulation  of  the  IRMC  in  the  azimuthal  mode 
domain  is  divided  into  two  parts.  In  the  first  step  the  modal 
vector  P is  used  to  calculate  the  actuator  signal  amplitudes 
| according  to  equation  (7).  For  brevity  the  corresponding 
transformation  is  represented  by  the  modal  transfer  matrix  T'4 
in  figure  4.  Equation  (6)  gives  the  secondary  modal  amplitudes 
S_m  generated  at  the  sensor  locations  with  help  of  the  modal 
transfer  matrix  f" . The  residual  amplitudes  Em  are  the  result 
of  the  superposition  of  the  primary  and  the  secondary  tone. 

The  presence  of  noise  leads  to  deviations  between  the  “meas- 
ured” and  the  actual  modal  amplitudes  at  each  sensor  location, 

eJxj)  = Pn(.xj)  + Smixj) 

« (8) 

h= 1 

For  that  reason  an  iterative  algorithm  is  used  in  the  following 
steps,  which  should  optimise  the  actuator  amplitudes  step  by 
step  in  order  to  minimise  the  residual  mode  amplitude. 

The  task  is  solved  by  a multiple-channel  LMS-algorithm  in  the 
azimuthal  mode  domain.  The  iterative  algorithm,  originally 


first  step 


iterative  steps 


Figure  4:  Initial  and  iterative  steps  implemented  in  the 
IRMC  algorithm  in  the  azimuthal  mode  domain. 
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designed  for  use  in  the  frequency  domain  (see  Kuo  and  Morgan 
[6]),  minimises  the  residual  mode  amplitude  according  to  equa- 
tion (3)  and  gives  the  update  prescription  of  the  modal  actuator 
vector  for  step  n as 

L(n  + = £,/«)  + U ■ (Ty'EJn),  (9) 

where  ji  is  the  adaptive  step  size  parameter. 

If  all  parameters  of  the  sound  field,  i.e.  frequency,  temperature 
and  mean  flow  velocity,  are  constant  in  time,  then  the  spatial 
structure  of  the  sound  waves  in  the  duct  repeats  periodically. 
Thus,  since  the  actuator  amplitudes  are  set  in  relationship  to  the 
trigger  of  the  rotor,  there  is  no  need  to  process  the  iterative 
algorithm  in  real  time. 

In  fact,  calculation  of  the  initial  modal  actuator  amplitudes 
| (0)  in  the  first  step  of  the  iterative  algorithm  is  not  manda- 
tory, however,  it  leads  to  a faster  convergence  of  the  iterative 
process. 

5.  Analysis  of  the  modal  transfer  matrix 

In  theory,  without  additional  noise,  the  modal  actuator  signal 
amplitudes  given  in  equation  (7)  should  produce  a cancelling 
tone  Sm  that  exactly  eliminates  the  primary  tone  P . However, 
this  is  true  only  for  azimuthal  modes  with  a low'  radial  content, 
which  is  equivalent  to  a corresponding  modal  transfer  matrix  of 
low'  rank.  In  this  chapter  it  will  he  shown,  that  with  increasing 
radial  mode  order  n numerical  restrictions  lead  to  limitations 
with  regard  to  the  solution,  w'hich  is  exact  analytically. 

All  calculations  presented  in  this  paper  have  been  carried  out  for 
the  configuration  depicted  in  figure  2.  A uniform  flow  ofM=0.3 
and  equally  spaced  rings  in  the  axial  direction  of  the  duct  arc 
assumed  for  various  parameters  k0ra  (k0  denotes  the 
wavenumber  for  free  field  propagation).  The  axial  distance 
between  two  consecutive  rings  is  denoted  by  & . 

Calculation  of  the  modal  transfer  matrix  Tj'  involves  the  inver- 
sion of  a hermitian  matrix,  named  U in  the  following.  The 
accuracy  of  the  modal  transformation  process  strongly  depends 
on  the  condition  of  U_  , whereby  a low  condition  corresponds 
to  a transformation  with  good  accuracy.  In  other  words,  the 
condition  of  U_  is  a measure  for  the  quality  of  the  detection  of 
the  modes  at  the  specified  positions.  Thus  by  an  analysis  of  71' 
resp.  U_  , the  optimum  positioning  of  the  sensors  and  actuators 
for  application  of  the  IRMC-technique  can  be  examined. 

In  figure  5 the  variation  of  the  condition  of  U_  with  ring  dis- 
tance Sx  is  shown  for  all  propagating  azimuthal  mode  orders 
and  for  three  different  values  k0ra  . Obviously,  the  lower  the 


mode  order  m the  smaller  the  intervals  of  Sx  with  a low  condi- 
tion are.  This  can  be  explained  by  the  fact,  that  the  lowest  mode 
order  m=0  has  the  largest  radial  content;  see  for  illustration 
figure  6,  where  the  normalised  axial  wavelengths  of  all  propa- 
gating modes  (m,n)  are  depicted  for  k0ra  -16.3.  For  a low'  mode 
order  m the  differences  between  the  corresponding  axial  wave- 
lengths /L^,  (n=0,...,N)  become  smaller,  which  makes  the  im- 
plicit modelling  of  the  radial  modes  by  the  mathematical  map- 
ping more  difficult.  Furthermore  the  complexity'  of  this  problem 
increases  with  increasing  k0ra  . On  the  other  hand,  axial  ring 


mode  order  m 


Figure  5:  Condition  of  the  hermitian  matrix  U„,  as  function 
of  the  axial  ring  spacing  &x  calculated  for  all  propagating 
mode  orders  m and  different  values  k0ra  . 
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spacing  is  almost  arbitrary  when  the  azimuthal  mode  order  m is 
large  and  the  radial  mode  content  low. 

In  general,  it  is  the  lower  mode  order  m that  impose  the  essential 
restrictions  with  regard  to  the  positioning  of  the  rings.  Compari- 
son of  the  examples,  shown  in  figure  5 leads  to  the  conclusion, 
that  for  simultaneous  control  at  several  frequencies,  the  proper 
choice  of  the  ring  positions  becomes  a critical  issue. 

6.  Analysis  of  the  IRMC-Algorithm 

The  performance  of  the  IRMC-algorithm  in  the  presence  of 
noise  is  examined  in  this  chapter.  The  calculations  have  been 
carried  out  for  the  same  parameters  as  in  chapter  5 for 
kBra  =16.3  . The  magnitude  of  the  amplitudes  of  the  primary 
modes,  that  are  radiated  by  the  rotor  and  propagate  upstream, 
are  randomly  distributed  about  = 85  dB.  The  reflected 
modes  propagating  in  flow  direction  are  assumed  to  be  in  the 
range  of  P*n  = 70  dB.  The  iterative  control  system  is  considered 
to  be  converged,  if  the  mean  residual  amplitude  over  all  sensors 
equals  50  dB. 


mode  order  m 


2.0  H — 1 


mode  order  m 


Figure  6:  Normalized  axial  wavelengths  Kmn  and  Xtm  of  the 
modes  (m,n)  propagating  in  resp.  against  flow  direction  for 
k0r  = 16.3. 


As  pointed  out  in  chapter  5,  a system  with  an  ill-conditioned 
transfer  matrix  is  very  sensitive  to  small  disturbances  of  the 
input  vector.  This  is  confirmed  by  figure  7,  which  shows  the 
number  of  iterations  until  convergence  of  the  iterative  algorithm 
depending  on  the  axial  ring  spacing  Sx  . The  calculations  are 
perfonned  for  mode  order  m=8  and  different  noise  levels,  given 
by  the  modal  signal  lo  noise  ratio  a . The  larger  the  distur- 
bances, the  larger  the  deviations  between  the  calculated  output 
of  the  modal  matrix  transform  and  the  optimal  solution,  as  it  is 
described  by  equation  (8).  In  the  absence  of  noise  or  with  very 
low  noise  levels  ( a >30dB),  the  IRMC  cancels  the  primary  tone 
completely  after  the  first  step  for  any  Sx  . For  er  <30dB  and  a 
ring  distance  Sx  corresponding  to  a quite  good  condition,  the 
control  is  stable  and  after  a few  iterations  the  algorithm  also 
converges.  But  with  increasing  noise  level,  and  ring  spacings 
corresponding  to  an  ill-conditioned  system,  the  system  becomes 
unstable.  In  these  cases  lots  of  iterations  are  necessary  until 
convergence  is  reached,  or  the  control  system  even  does  not 
converge  at  all. 

In  Figure  8 a control  attempt  for  the  mode  order  m=8  and  an 
optimal  ring  spacing  Sx  is  shown  for  a - 5dB  in  comparison 
with  a A30dB.  In  the  coordinate  system  used  here,  the  rotor  is 
located  at  x=0  and  the  uniform  flow  goes  in  the  positive  axial 
direction.  Depicted  are  the  magnitudes  of  the  complex  azimuthal 
mode  amplitudes  over  the  duct  axis  after  the  iterative  algorithm 
has  converged.  Obviously  the  actuators  reproduce  the  primary 
mode  well,  resulting  in  the  cancellation  of  P_  . The  residual 
mode  amplitude  Em  strongly  depends  on  the  uncertainty  of  the 
measured  azimuthal  mode  amplitudes. 

The  contour  plot  in  figure  9 shows  the  residual  modal  ampli- 


Figure  7:  Number  of  iterations  until  convergence  as  a func- 
tion of  the  sensor  ring  spacing  for  various  modal  signal  to 
noise  ratios  cr  (in -8,  n=0,l,  k0ra  = 16.3). 
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tudes  Em  (x,  ra ) along  the  duct  axis  for  the  mode  order  m=8 
and  <7  =15  dB  for  different  ring  spacing  Sx  . As  examination  of 
the  matrix  U_  predicted,  for  a ring  spacing  Sx  connected  with 
bad  condition  of  U_  , the  system  does  not  cancel  the  primary 
mode  or  even  diverges,  which  is  indicated  by  residual  modal 
amplitude  Em  larger  than  the  modal  amplitude  Pm  of  the  pri- 
mary tone. 

Note,  that  control  of  the  primary  mode  at  the  locations  of  the 
sensors  is  successful  for  arbitrary  ring  spacings  Sr  : the  reduc- 
tion at  the  axial  positions  of  the  sensor  rings  are  clearly  visible 
in  the  contour  plot  as  white  lines  (the  arrows  mark  the  moving 
positions  of  the  rings  with  increasing  <&  ).  The  important  con- 
clusion is,  that  a control  observed  at  the  sensors  does  not  neces- 
sarily mean  a reduction  of  the  primary  mode  in  the  entire  duct. 


Figure  8:  Azimuthal  mode  amplitude  of  the  primary  and  the 
residual  tone  along  the  duct  axis  for  modal  signal  to  noise 
ratios  of  a =5  and  30  dB  after  convergence  of  the  fRMC 
(m=8,  n * =0,1)  Ay;,  = 16.3,  rotor  at  x=0). 


Figure  9:  Residual  azimuthal  mode  amplitude  along  the  duct 
axis  with  a modal  signal  to  noise  ratio  of  a =15  dB;  the 
arrows  mark  the  positions  of  the  sensor  rings  depending  on 
Sx  (m=8,  n =0,1,  k0ra  = 16.3,  rotor  at  x=0). 


The  upper  graph  of  figure  1 0 shows  a control  attempt  for  mode 
order  m=l  and  a =25dB.  The  azimuthal  mode  contains  10 
radial  modes  (cp.  figure  6)  and  should  be  detected  by  10  sensor 
rings  with  an  axial  spacing  Sx/ra  = 0.1  corresponding  to  a 
matrix  U_  with  minimum  condition.  Nevertheless  the  condition 
for  this  configuration  is  quite  bad  and  convergence  is  not 
reached.  This  example  again  confirms,  that  control  achieved  at 
the  sensor  location  does  not  imply  control  at  all  duct  locations. 
The  situation  is  improved  with  additional  sensor  inputs,  e.g. 
with  13  sensors,  which  is  shown  in  the  bottom  graph  of  figure 
10.  The  corresponding  matrix  U_  still  is  bad  conditioned. 
However,  if  the  inverting  process  of  the  hermitian  matrix  is 
carried  out  via  a singular  value  decomposition,  columns  corre- 
sponding to  singular  values  close  to  zero  can  be  rejected,  which 
results  in  an  improved  condition  and  performance  of  the  IRMC. 

Since  the  IRMC-technique  makes  use  of  wall  flushed  sensors 
and  actuators,  one  has  to  prove  that  not  only  the  modal  ampli- 
tudes at  the  duct  wall  are  controlled  but  also  the  whole  azi- 


Figure  10:  Azimuthal  mode  amplitude  of  the  primary  and 
the  residual  tone  along  the  duct  axis  for  a modal  signal  to 
noise  ratio  of  25  dB  and  10  sensor  rings  (top)  resp.  13  sensor 
rings  (bottom)  (axial  ring  spacing  <5x/'r  = 0.1,  m=l, 
n * =0,...,4,  kur:  = 16.3,  rotor  at  x=0). 


19-7 


muthal  and  radial  structure  across  the  duct.  Evidence  to  that 
regard  is  given  in  figure  1 1 and  12,  which  show  the  pattern  of 
the  primary  and  residual  azimuthal  mode,  P and  Em  resp., 
across  the  duct  corresponding  to  the  control  attempt  shown  in 
figure  10. 

7.  Summary 

The  results  of  the  analysis  and  simulation  show,  that  the  IRMC- 
concept  enables  an  effective  control  of  the  azimuthal  and  radial 
structure  of  acoustical  modes  in  circular  and  annular  ducts  with- 
out the  need  to  determine  the  radial  amplitudes  directly.  All 
azimuthal  modes  of  a tone  can  be  reduced  by  using  the  same 
ring  spacing  for  the  sensor  rings  and  actuator  rings.  In  this  con- 
text, the  lower  azimuthal  mode  orders  m impose  the  most  strin- 
gent requirements  for  the  axial  spacing.  For  effective  control  of 
an  azimuthal  mode  m with  a high  radial  content  (N  radial 
modes),  more  than  N sensor  rings  are  necessary.  Control  at- 


azimuthal  position  /"/ 

Figure  11:  Pattern  of  the  azimuthal  mode  amplitude  (mode 
order  m=l)  of  the  primary  tone  in  the  cross  section  of  the 
duct  at  x/ra  =-4.5  (parameters  as  in  figure  10). 
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tempts  at  several  frequencies  simultaneously  require  a detailed 
analysis  with  respect  to  the  optimal  ring  spacing. 

An  important  result  with  regard  to  the  experimental  analysis  of 
active  control  concepts  is,  that  noise  control  achieved  at  the 
positions  of  the  sensors  does  not  necessarily  mean  a complete 
control  at  all  positions  in  the  duct. 

The  simulation  results  show  further,  that  it  is  not  necessary'  to 
place  the  sensor  rings  in  the  axial  direction  according  to  a spatial 
Nyquist  criterion.  This  can  be  explained  by  the  fact,  that  in  the 
adaptation  process  only  the  radial  mode  amplitudes  belonging  to 
certain  discrete  modal  wavenumbers  have  to  be  fitted  im- 
plicitly. 

Another  result  is,  that  complete  control  is  possible  in  the  duct 
although  the  sensor  array  covers  less  than  half  of  the  maximum 
modal  wavelength. 

In  future  investigations  non  equidistant  axial  spacing  of  the  rings 
has  to  be  examined  with  respect  to  matrices  with  improved 
conditions.  Additionally  improved  methods  for  matrix  inver- 
sions have  to  be  tested.  A perturbation  method  may  give  more 
information  in  this  context. 
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Figure  12:  Pattern  of  the  azimuthal  mode  amplitude  (mode 
order  m=l)  of  the  residual  tone  in  the  cross  section  of  the 
duct  at  x/ra  —4.5  (parameters  as  in  figure  10). 
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